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corrections are exponentially small in the volume, but can become pronounced for small twist 
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magnetic moment and the electromagnetic radii from simulations at zero lattice momentum. Our 
formalism, moreover, can be applied to any nucleon matrix elements. 

PACS numbers: 12.38.Gc, 12.39.Fe 



* fj j iang@itp . unibe . ch 
^ bctiburz@umd.edu 



1 



I. INTRODUCTION 



Understanding QCD in the strongly interacting regime remains a challenging problem in 
physics. Simulations of QCD on Euclidean spacetime lattices are making progress towards 
a quantitative understanding of the non-perturbative dynamics in QCD [l|. Lattice QCD 
simulations usually employ periodic boundary conditions for the quark and gluon fields. 
Consequently the available hadron momenta are limited to periodic momentum modes of the 
lattice, k = 2irn/L, where n is a triplet of integers and L is the lattice size in each of the three 
spatial directions. On typical lattices, the smallest available lattice momentum is about 400— 
500 MeV. This presents a severe limitation for the study of observables appearing in matrix 
elements at low momentum, and low momentum-transfer. At present, such observables 
cannot be investigated directly using periodic boundary conditions, and models are used to 
perform momentum extrapolations. 

For large enough volume, the physics should be independent of the choice of boundary 
conditions. There is freedom in choosing boundary conditions for fields; however, the action 
must be single valued so that observables are well-defined. For a generic matter field 0, we 
can impose a twisted boundary condition in the z-th direction of the form, see e.g. 0], 

4>(xi + L) = U 4>(xi), 

where U is a symmetry of the action and XPU = 1. For the quark flavors in QCD, the 
diagonal flavor rotations can be used to implement what are called flavor-twisted boundary 
conditions. With U of the form U = exp(i6i), the matter field <ft has kinematic momentum 
k = (2-Kn + 0)/L which can be varied continuously by choosing different values for 0. 
The ability to produce continuous hadron momentum has made flavor-twisted boundary 
conditions attractive to lattice QCD @, |, H, 0, 0, 0, M, M, Q G3, M, H H EE 0, M ■ 



In this work, we detail the finite volume modifications to nucleon form factors of the vector 
current. We use twisted boundary conditions on the active quarks in the current insertion, 
and, of course, are limited to only connected contributions from the current. Heavy baryon 
chiral perturbation theory is utilized to estimate the volume dependence of nucleon current 
matrix elements. Let us summarize our main findings. 

• Finite volume modifications can be sizable especially for the magnetic contribution, 
and for small twist angles. 

• The use of Breit frame kinematics does not dramatically reduce or simplify the finite 
volume corrections. The volume effect for magnetic observables in the Breit frame 
roughly doubles compared to the rest frame. This situation is unlike the meson vector 
current tlj 



With twisted boundary conditions on only the active quarks, the finite volume correc- 
tions depend on an unphysical and unknown parameter, g 1 . This dependence arises 
as an artifact of the enlarged valence flavor group, and a lattice determination of g\ 
would help in accounting for volume corrections. 

Results obtained here are qualitatively similar to those obtained from isospin-twisted 
boundary conditions for the nucleon isovector form factors [13]. In that method, va- 
lence li-quarks are twisted differently than the valence ci-quarks without introducing 
extra fictitious flavors. We show, moreover, the flavor symmetry employed by that 
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method can be used to eliminate the dependence on g\. Consequently finite volume 
corrections can be reliably estimated for that case in terms of known low-energy con- 
stants. 

Our presentation has the following organization. First in Section [Til we detail the flavor 
twisted boundary conditions, and incorporate them into heavy baryon chiral perturbation 
theory. We show how the graded group SU(7\5) accommodates twisting of the active va- 
lence quarks in the baryon sector. In Section [TTTJ we compute finite volume corrections to the 
nucleon mass, and derive the induced mass splittings due to flavor twisted boundary condi- 
tions. Numerically these splittings are estimated to be at the percent level or less on current 
lattices. Finite volume corrections to the vector form factors of the proton and neutron are 
determined in Section HVl (complete expressions are given in Appendix |AJ. These results ap- 
ply to the connected contributions allowing access to isovector quantities, but not isoscalar. 
We show that terms arising from broken cubic invariance can lead to non-negligible volume 
effects in the region of small twist angles. Results for rest frame and Breit frame kinematics 
are compared. The Breit frame does not offer any substantial advantages with respect to 
volume effects. Complete results for isovector current matrix elements at finite volume using 
the method of isospin twisted boundary conditions are displayed in Appendix [B] For rest 
frame kinematics, these results are shown to be independent of the unphysical parameter 
g\. Appendix O collects functions and identities useful for the evaluation of finite volume 
effects. Finally a brief summary concludes our work. 

II. FLAVOR TWISTED BOUNDARY CONDITIONS AND BARYON CHIRAL 
PERTURBATION THEORY 

To address the consequences of twisted boundary conditions in lattice calculations of 
baryon properties, we describe the underlying effective theory in the baryon sector. First 
we detail the partially twisted boundary conditions employed. Next we include these effects 
in chiral perturbation theory, and then heavy baryon chiral perturbation theory. 

The quark part of the partially quenched QCD Lagrangian is given by 

12 

C= Q j (P+ m Q)j k Qk- (i) 

j,k=l 

The twelve quark fields transform in the fundamental representation of the graded group 

«S77(7|5), and appear in the vector Q T = (u ,ui,u 2 ,di,d2,jJ,Uo,Ui,U2,di,d2)- In addition 
to the valence uq, u\, «2, d\, and d 2 quarks, we have added their ghost quark counterparts uo, 

Ui, u 2 , d~i and d 2 , which cancel the closed valence loops, and two sea quarks j and I. In the 
isospin limit, the quark mass matrix of SU(7\5) reads ttiq = diag(m u l 5x5 , rrijl 2x2 , m M l 5x5 ) 
in block diagonal form, where the blocks correspond to valence, sea and ghost sectors. QCD 
quantities can be recovered in the limit rrij — > m u . The additional up and down quarks 
are fictitious flavors differing only by their boundary conditions. There is one more up-type 
quark than down-type quark because we focus on a theory that will yield proton matrix 
elements. Neutron matrix elements can always be derived trivially by interchanging up and 
down charges in the final result. 1 



To consider both proton and neutron properties in the same theory, we would need to enlarge the flavor 
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The hats denote fields satisfying twisted boundary conditions. We require that the quark 
fields satisfy boundary conditions of the form 

Q(x + Le r ) = exp (i9^T a ) Q(x), (2) 

where e r is a unit vector in the r th spatial direction, L is the spatial size of the lattice, and 
the block diagonal form of the supermatrices T a is 

T a = diag(T a ,0,T a ). (3) 

Here we choose T a to be generators of the U(5) Cartan subalgebra. Notice that by Eq. ([3]), 
the sea quarks remain periodic at the boundary. This reflects a partially twisted scenario. 
Twist angles can be changed without necessitating the generation of new gauge configu- 
rations, because the fermionic determinant, which arises solely from the sea sector, is not 
affected by the twisting. 

Redefining the quark fields as Q T = (uo, Ux, u 2 , d%, d 2 , j, I, uq, «i, u 2 , d\, d 2 ), with Q(x) = 
V jf (x)Q(x), where V(x) = exp(i6 a ■ xT a /L), we can write the partially quenched QCD 
Lagrangian as 

12 

c= J2Q j (i>+ m o) k Qk, (4) 
j,k=l 3 

where all Q fields satisfy periodic boundary conditions, and the effect of twisting has the form 
of a uniform gauge field: = D^ + iB^, where = (6 a T a /L, 0). It will be easier to treat 
the twisting in the flavor basis of the valence and ghost sectors rather than in the generator 
basis, thus we write a T a = diag(0, 8 U , /u , d , 0' d ), and similarly for B^ which appears 
as 5 M = diag( J B™ 1 ,0,5™ 1 ) in block diagonal form, with Bjf = diag(0, B%, B'f, B^ B'fl. 
Momentum transfer will be generated using flavor changing currents from u\ to u 2 , or from 
d\ to d 2 . Notice we keep the uq quark periodic; it will play the role of spectator. 

The low-energy effective theory of QCD is chiral perturbation theory, which describes 
the dynamics of pseudoscalar mesons emerging from spontaneous chiral symmetry breaking. 
The mesons of partially quenched chiral perturbation theory pji, 20, 2l|, 22, 23] are described 



by the coset field E, which satisfies twisted boundary conditions. This field can be traded 
in for E, defined by E(x) = V* '(x)E(x)V (x) , which is periodic at the boundary j^]. In terms 
of this field, the Lagrangian of partially quenched chiral perturbation theory appears as 



C = f —sti [p^D^ - X str (m^E + E^q) + (5) 

The action of the covariant derivative is specified by -D^E = -D^E + i[i? /t ,E]. The 
parameter / is the chiral limit value of the pion decay constant, and in our normalization, 
/ = 0.13GeV. The above Lagrangian contains only the lowest-order terms in an expansion 
in quark mass rriQ, and meson momentum k 2 . The periodic meson fields contained in the 
twelve-by- twelve matrix are realized nonlinearly, E = exp(2i</> / '/) . The matrix <fi has the 
form 

( M VV M vs xL\ 

M sv M ss X \s ■ (6) 

Xgv Xgs MggJ 



group further to SU(8\6). 
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The mesons of M vv (M gg ) are bosonic and are formed from a valence (ghost) quark-antiquark 
pair. These matrices have the form 



{ Voo Vol Vo2 K 01 ^(^A 

VlO Vll Vl2 *li 7Tl2 

^20 ^21 ??22 7r 21 ^22 

^[o *Ti n u Vn Vf 2 

\ n 20 n 21 n 22 V21 V22/ 



, and M gg 



\ 



VoO Vol V02 71 01 7r 02 

VlO Vll Vl2 111 ^l + 2 

V20 V21 V22 ^21 n 22 

kJx frfa V11 V12 

\ 7T 20 n 21 n 22 V21 V22) 



The r\\- (77^) mesons have quark content 77. 

have quark content tt^ ~ U{dj (fc^ 
and contained in M vs (M ss ) as 



' ~ QiQj (vli ~ <M 7 )> while the 7r^ (7r+) mesons 



VI13 111] /I — 

■UjC?, ). The valence-sea (sea-sea) mesons are bosonic 



M x , 



(fijuo 'rjui 4 > ju2 tfijdi 4 > jd2 
4>lu <f>lui 4>lu 2 4>ldi </>Zd 2 



and Af, 



Vj Kji 
7Ty 77; 



Mesons contained in (Xgs) are built from ghost quark, valence antiquark (sea antiquark) 
pairs and are thus fermionic. These states appear as 







feoui 


fe U2 


feodi 


fe d 2 ^ 






fe i\ 




0mi«o 






feldl 


4 > uid2 




< fiuij 


fail 


Xgv 




^rU 2 Ui 


4 > U2U2 


Vuzdi 


4 > U2d2 


, and Xgs = 




fe 2 l 




feiu 


feim 


felU 2 


feldl 


feid 2 




<t>dlj 


fell 






Vd2U X 


Vd 2 U2 


§d,2dx 


feldj 




\<t>dv 


fell] 



Expanding the Lagrangian in Eq. 
content QQ' have mass 



to lowest order, one finds that mesons with quark 

(7) 
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2 ( m Q + m Q' 



Thus in infinite volume all mesons fall into one of three groups of mass degenerate states: 
valence-valence pions m 2 = 8\m u /f 2 , valence-sea mesons m 2 u = 4A(m u + rrij)/f 2 , and 
sea-sea pions m|- = 8\rrij/ f 2 . In partially quenched simulations, one measures the valence- 
valence and sea-sea pion masses. The valence-sea mass is given by the average of the other 
two, up to possible discretization errors that arise in hybrid actions. 

The flavor singlet field, $0 — ^^str 0, additionally acquires a mass /Jo which arises as a 
consequence of the U(1)a anomaly. Taking this mass to be large, the flavor singlet field is 
then integrated out of partially quenched chiral perturbation theory, however, the propaga- 
tors of the flavor- neutral fields deviate from simple pole forms 22|, |23j . There are two useful 
simplifications to note: twisted boundary conditions have no effect on the flavor neutral sec- 
tor, and all valence-valence flavor neutral states are degenerate with mass m n . For a,b = uq, 
ui, U2, di, or G? 2 , the leading-order r] a r]b propagator is thus given by 



GvaVb = S < 



k 2 + m 



ab 



k 2 + m 2 2 (k 2 + m; 



2^2' 



The flavor neutral propagator can be conveniently rewritten as 

Sab - . ( I 



>VaVb 



k 2 + m 2 



+ n 



ab 



k 2 + m 2 



(8) 



(9) 
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where 



Hat (A) 



1 

2 



l + (mi 



m 



33' 



d 
dm? 



A. 



(10) 



To include baryons into partially quenched chiral perturbation theory, one uses rank 
three flavor tensors @, HI, 0, HJ. In SU(7\5), the spin-i baryons are described by the 

572-dimensional supermultiplet B l ^ k , while the spin-| baryons are described by the 300- 
dimensional supermultiplet T^ k (28[ . The baryon flavor tensors are twisted at the boundary 
of the lattice. In the r th spatial direction, both tensors satisfy boundary conditions of the 
form [7J] 



B ijk {x + e r L) 



set r a 



JJBZ T a 
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kk 



Bijk(x) 



Thus we define new tensors B^ and T^ k both having the form 

B ijh (x) = Vi(x)vUx)Vl(x)B ijk (x). 



(11) 



(12) 



These baryon fields satisfy periodic boundary conditions, and their free Lagrangian has the 
form 



C = -i [ Bv^V^B 



2a { M Q) (BBM+) - 2[3<£ Q) {BM+B) - 2a ( ^ Q) (BB) str (M 



i(PQ) 



(PQ) 



{T^VfX) + A (T v %) + 2^ Q) (T u M + %) + 2a^ Q) (7 v %) str (M 



(13) 



where = (0, 0, 0, i) is the Euclidean four-velocity in the rest frame. The mass operator 
M. + is defined by A4 + = | (^rriQ^ + £?^q£), with £ = yS, and the covariant derivative 
acts on B and 7^ fields in the same manner, namely 



ijk. 



X 



l (B; + Bj 1 + B k )B^ k [ 



x 



(14) 



with 



[£> M £p = d^B ijk + {V^fB l3k + (-)^ 3 +m)(y^yi B uk + ^_Y^+v 3 )(v k +m)^) kl B ijl , (15) 

where the vector field of mesons is given by = ~ (^D^ + ^D^^j. This free La- 
grangian contains a number of low-energy constants but has precisely the same form as in 
the SU(4\2) partially quenched theory. Restricting the baryon multiplets to the sea sec- 
tor, so that all flavor indices are either 6 or 7, we have nucleons and deltas made only of 
sea quarks. Hence the matching conditions are precisely the same as those used to match 
SU(4\2) onto SU(2) by restricting the former to the sea sector. The relations between the 
low-energy constants appearing in Eq. (fT3|) and those of SU(2) chiral perturbation theory 
will not be needed here, but are given in [271 ] . 

The leading order partially quenched interaction Lagrangian between the baryons and 
mesons appears as 



C= 2n ( BS^BA^j + 2(3 (BS^B 



2H(T u S fl A fl %)+\l-C 



r v A v B 



BA V % 



(16) 
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where the effects of partial twisting show up in the axial-vector field of mesons = 

| (^D^ — ^Dp^j. The interaction Lagrangian has the same form as the S77(4|2) theory 

of baryons, hence the matching conditions to 577(2) are identical. The familiar low-energy 
constants of SU(2) are identified as follows j27|: — |a — gAN = — C, and #aa = 7~L. 
Notice there is an extra free parameter in the partially quenched interaction Lagrangian 
compared to that of SU(2) chiral peturbation theory. We shall write our expressions in 
terms of g^ and the combination g\ — \a-\- Dependence on gi must drop out in the 
QCD limit, which, for the case at hand, requires both mj — > m u and L — > oo. 



III. NUCLEON MASS 

To begin, we determine the nucleon mass in the presence of partially twisted boun dary 



conditions. In the isospin limit of S77(4|2), the proton and neutron are degenerate [271. 129 
Flavor twisted boundary conditions, however, break the valence flavor symmetry, hence the 
nucleons are no longer degenerate. The nucleon mass splittings arise from finite volume 
effects induced by the boundary conditions. Effects of this type can be treated using chiral 
perturbation theory at finite volume. We work in the p-regime throughout, where m n L 3> 1 



so that zero modes of the pion field do not become strongly coupled [3Cj, l31|, [32] . We estimate 
the size of the mass splittings on current-sized lattices, and show that their effect can be 
neglected for the determination of nucleon observables using twisted boundary conditions. 

In the infinite volume limit with held fixed, the nucleon mass is unaffected by the 
boundary conditions. This follows from a generalization of the argument presented for 
mesons in jBf. One merely realizes that the nucleon propagators are not boosted in heavy 
baryon chiral perturbation theory, because v^B^ = 0. The remaining momenta in a given 
diagram are mesonic, and are boosted according to flavor. Now since there are no flavor 
changing interactions, the sum of boosts at each vertex is zero. This, along with v^B^ = 0, 
assures us we can always shift loop momenta to cast any diagram into a form where the only 
^-dependence is that from external momenta. These contributions should be thought of as 
kinematical rather than effects which arise in the loops from chiral dynamics. The mass does 
not receive dynamical corrections from the boundary conditions, but the energy depends on 
the external momentum and has a kinematic dependence on the boundary conditions, e.g. for 
a nucleon with lattice momentum k 

^ = M » +1 W + -- (17) 

where B = 0/L, and ... denotes terms that are higher order in 1/Mjy. Here only one 
valence quark in the nucleon has been twisted, and by an angle 6. When we try to apply the 
same reasoning at finite volume, no shifts of the internal momenta are possible because the 
loop momenta are discrete, while the twisting parameters are continuous. Thus there is a 
dynamical dependence on the twisting parameters arising from loops, and in finite volume, 
M N will depend upon B. 

The mass of the nucleon in the chiral expansion can be written in the form 

M N = MM - M$(ja) - M^ 3/2) (/i) + . . . , (18) 

where /x is the renormalization scale, and Mjy denotes the contribution to the nucleon mass 
of order m" . The linear quark mass dependence arises from the local operators in Eq. ffT3|) 
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FIG. 1: Diagrams contributing to the nucleon mass and wavefunction renormalization in partially 
quenched chiral perturbation theory. A thin (thick) line denotes a spin- 1/2 (spin-3/2) baryon, 
while a dashed line denotes a meson. Partially quenched hairpins are depicted by a crossed dashed 
line. 



at tree level, while the leading non-analytic contribution 0{rr?J 2 ) arises from the one-loop 
diagrams shown in Figured] The local interactions do not contribute to finite volume effects, 
only the meson loops that are shown in the figure. For periodic boundary conditions, the 
finite volume effects on the nucleon mass have been determined in (33[. To express the finite 
volume corrections to the nucleon mass with flavor twisted boundary conditions, we require 
the mode sum 



JC(m, B, A) 



-y 

7-3 



(k + B) 



L 3Z^ [{k + B) 2 +(3 2 A 



2 13/2 



dk 



k 2 



(2vr) 3 [fc 2 + /? 2 ] 3 / 2 



(19) 



with (3\ = A 2 + 2AA + m 2 , and k = l-KnjL where n is a triplet of integers. Evaluation of 
this function, as well as other finite volume sums, is discussed in Appendix S 

Consider purely valence nucleon states with exactly one twisted quark. These will be 
the only nucleons relevant in the computation of matrix elements with twisted boundary 
conditions. It is easiest to classify these states according to their representations under the 
valence subgroup of the two degenerate untwisted quarks. In our formulation, we must set 
the twist angles for these quarks to zero by hand. There is both a singlet, 1, and triplet, 
3, representation for singly twisted nucleons under the untwisted valence SU{2). For the 
mass of a nucleon with one twisted valence quark in the 3 representation, we find the finite 
volume shift SM Na (B) is given by 



SM N JB) 



9nN 3 NM m ^ °> °) + g'nNsNs^^n, B, 0) + g juNaNa K,{m jui 0, 0) 



2 J2 j ^irNsN3 

+ g'jlN 3 N 3 fc( m ju, B, 0) + (g A + gi) 2 H uu (jC(m n , 0, 0) 

1 



+ q9an 



K{m v , 0, A) + 5!C(m n , B, A) + 2)C(m ju , 0, A) + AJC(m ju , B, A) 



(20) 



S 



The effective axial couplings are defined by 



9ttN 3 N 3 



a' 2 



9juN 3 N 3 
12 

9juN 3 N 3 



-{g 2 A + 2g A g 1 +g 2 l /A), 



{g\- 9a9i -5#i/4), 
(±g A + 2g A g 1 + g 2 1 ), 



1 2 

2 91 ' 



(21) 



On the other hand, for a nucleon in the 1 representation, we find the finite volume shift 
SM Nl (B) given by 



SM Nl (B) 



^jViiVx^^Tr, 0, 0) + g'* mNl K,{m^ B, 0) + g% NlNl )C(m ju , 0, 0) 
+ 9jlNiN 1 J ^( m ju, B, 0) + (g A + gi) 2 H uu \K.(m n , 0, 0) 
JC(m n , 0, A) + IC(m n , B, A) + 2}C(m ju , 0, A) 



+ ^9 AN 



The effective axial couplings for the singlet nucleon mass are defined by 



(22) 



9l NlNl = 7:(9A-^9A9i-ngl/4=), 



a' 2 

9 ttNi N± 



9 

-(5^ + 7^-^/4), 



9]uN lNl = ^9 2 A + ^9A9i + 7gl), 



n' 2 



(8g 2 A + 4g A g 1 + g 2 1 /2) 



(23) 



In the limit B 0. there is no difference between the representations, and we recover 
accordingly the finite volume shift of the partially quenched nucleon mass [33]. In the rest 



of this section, we work for simplicity at the unitary mass point rrij = m u , so that m 2 u 

We consider nucleon splittings for two cases that are of interest in current matrix elements: 
rest frame kinematics, and Breit frame kinematics. In the rest frame kinematics, the initial 
nucleon is at rest, and hence completely untwisted. The final nucleon has been given a 
boost by twisting one of the quarks by 6. In this case, there are three mass splittings among 
the various nucleons: that between the 3 and untwisted nucleon, that between the 1 and 
untwisted nucleon, and that between the 3 and 1 nucleons. The relative change in these 
splittings is given by 



AMq 



M N3 (B)-M N 



M, 



N 



(24) 



2f 2 M, 



N 



+ 9 AN 



{9 a ~ 9a9\ + £i/4) lC(m n , B, 0) - JC(m n , 0, 0) 



JC{m n ,B,A) - JC(m n , 0, A) 



(25) 
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A Mi 



M N JB)-M, 



N 




(26) 



(13^ + llg A9l + g 2 jA) K{m^ B, 0) - JC(m n , 0, 0) 



K,(m 7T , B, A) — JC(m n , 0, A) 



(27) 



and 



AM 3 _i 



M N AB)-M N AB) 




(28) 



(10^ + Ug A9l - gl/2) K{m„, B, 0) - JC(m % , 0, 0) 



lC(m n} B,A) - /C(m^,0, A) 



(29) 



respectively. 

On the other hand, for the Breit frame kinematics the initial state nucleon has one 
quark twisted by 6, while the final state nucleon has one quark twisted by —6. The finite 
volume modification given by the function /C(m, B, A) in Eq. (Tl9|) is even with respect to B. 
Thus initial and final state nucleons in the same representation of the untwisted SU(2) are 
degenerate. The only non-vanishing splitting is between the different representations, but 
on account of evenness in B, this splitting is identical to AM 3 _! given above in Eq. (1281) . 

Numerically we can estimate the nucleon splittings by using phenomenological input for 
the low-energy constants: g A = 1.25, gAN = 1-5, A = 0.29 GeV, Mn = 0.94 GeV, and 
/ = 0.13 GeV. For the unknown partially quenched axial coupling g±, we assume it is of 
natural size and vary it within the rage —2 < g\ < 2. Each of the mass splittings AM3, 
AMi, and AM^-i is a maximum when = 7r(l, 1,1). We choose this value for 6 to 
investigate the worst case scenario. In Figure we investigate the relative mass splittings' 
dependence on L for a fixed value of m n , which is chosen to be 0.25 GeV. Values shown 
for the maximal splittings are all less than five percent. We will thus neglect the nucleon 
splittings in our analysis below. 2 



IV. NUCLEON ISOVECTOR FORM FACTORS 



Electromagnetic form factors appear in vector current matrix elements of the nucleon. 
In terms of Dirac and Pauli form factors denoted by F^(Q 2 ) and FJ?(Q 2 ), respectively, the 
nucleon current matrix element has the decomposition 



(N(P') 



\J™\N(P)) = u(P') 



2M, 



N 



U(P), 



(30) 



2 Additionally partially twisted isospin splittings in the meson sector have been shown to be negligible on 
current sized lattices 



13[. The same is true of the infrared renormalization of the twist angles. These 



effects will hence also be neglected. 



10 



m n = 0.25 [GeV] 



m„ = 0.25 [GeV] 




L [fm] 

FIG. 2: Numerical estimates for the maximal nucleon splittings at finite volume, AM3, AMx, 
and AA^3_i given in Eqs. (pHJ) . ([26]) . and ([28]) . respectively. We plot each relative splitting as a 
function of L with the lattice pion mass fixed at = 0.25 GeV. The twist angles are fixed at 
= 7r(l, 1, 1) to give the maximal splittings. The bands arise from variation of the parameter gi 
assuming naturalness. 



where = {P' — P) M is the momentum transfer. In QCD, the electromagnetic current is 
given as J^ m = q u u^^u + d^^d. In heavy baryon chiral perturbation theory, the decom- 
position of the current appears in terms of the Sachs electric and magnetic form factors, 
G^(Q 2 ) and G^(Q 2 ), i.e. one has 



(N V (P')\J™\N V (P)) -- 
with the relations 

G N E {Q 2 ) 
G N M {Q 2 ) 



u v 



v^G N E {Q A ) - 



M, 



G N M {Q 2 ) 



N 



U, 



F N (Q' 2 

N 1 



F\{Q) + F"(Q 



(31) 

(32) 
(33) 



We have appended velocity subscripts in Eq. ( 1311 for clarity. The u v are two-component 
Pauli spinors. 

To calculate these form factors with twisted boundary conditions on the lattice, one 
writes the current matrix element in terms of the various quark contractions with the elec- 
tromagnetic current. The propagators coupling to the current in each contraction we call 
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the active quark propagators. These are the propagators determined with twisted boundary 
conditions. Omitted from this calculation are the current insertions on quark lines that are 
self-contracted. These disconnected contributions are notoriously difficult to calculate us- 
ing lattice QCD. This difficulty notwithstanding, their contributions cannot be modified to 
produce continuous momentum transfer between the initial and final state hadron. As with 
present-day lattice calculations, we too will omit these contributions but with the caveat 
that their eventual inclusion will be limited to hadrons with Fourier momentum modes of 
the lattice. 

Now we discuss precisely how to calculate the connected part of the nucleon form factors 
in the effective theory. To specialize to the application of twisted boundary conditions on 
the active quarks, we must separate the current into two pieces, 

Jl = QuUiluV^ (34) 

(35) 

By evaluating matrix elements of J* with U = 6, ,u = 0', and J 2 with d = 6, ,d = 6', 
both currents induce momentum transfer from P = 6/L to P' = O'/L. First let us consider 
proton matrix elements. Considering the quark-level contractions, we find 



-{N 1 (uu 2 d 1 



e d =o 



+ -j\Na{uu 2 di] 



\ J l\ N * 



uuidx)) 



e d =o 



+ {N^UUd 2 )\jl\N^UUd X )) ^ (p(P')\J™\p(P)}conn C ct cd . (36) 

The subscripts on N refer to the representation under untwisted isospin, and parenthetically 
we list the quark content. We treat the active quark twists as implicit: each is from an initial 
quark, u\ or d\, with twist 6 to a final quark, w 2 or d 2 , with twist 6' . We stress that Eq. fl36|) 
provides the recipe for calculating in the effective theory what is implemented on the lattice 
by twisting the active quarks. For the neutron there is a similar construction, however, it 
is easiest to appeal to charge symmetry (isospin rotation by vr/2) from which follows the 
relation 



(n(P')\J™\n(P)) = (p(P>)\J™\p(P)) 



(37) 



In partially quenched QCD, the current is defined by J° = Q^^Q a Q. The choice of 

supermatrices Q a used to extend the charges is not unique [34]]. One should choose a form 
of the supermatrices that maintains the cancellation of valence and ghost quark loops with 
an operator insertion (35|. For the flavor changing currents we consider, the simplest choice 



Q l )ij = qu5i 3 5 j2 and (Q 2 



Qd 8i5 8j4 results in the correct physics. This is because 



operator self-contractions automatically vanish; thus, any non-zero charges in the ghost 
sector must ultimately yield zero, and consequently would be superfluous. Charges in the 
sea, while not superfluous, are absent due to restricting to the connected part of three-point 
functions. 

Operators that contribute at tree-level to the electromagnetic currents in SU(2) chiral 
perturbation theory are contained in the Lagrangian 



C 



J - (N[S„ S V ]N) tr(Q)i^ - J|L (N[S„ S V ]QN) F f 



2M 

~(NN) tr(Q)v,d u F^ 

x 



fiv 



£ (NQN) v,d u F, u , 
x 



(38) 
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where Q = diag(g„, qa) is the electric charge matrix. The operators with coefficients liq 
and fx i give the leading local contributions to the magnetic moments, while the operators 
with coefficients cq and cj give the leading local contributions to the electric charge radii. 
Operators for the magnetic radii occur at one higher order than the leading loop contri- 
butions. The combination |/i + is isoscalar, while /xj is isovector. Analogous linear 
combinations of cq and cj form isoscalar and isovector contributions to the charge radius. 
In writing down the analogous terms in the partially quenched chiral Lagrangian, Q — > Q a . 
Because of the condition str(Q a ) = 0, we see that there will be missing information in the 
partially quenched theory: there will only be operators where Q a transforms under the ad- 
joint, because of the lack of singlet component. This is the effective theory manifestation 
of neglecting the disconnected contributions. Consequently we will not be sensitive to the 
isoscalar combination low-energy constants. 

To consider the baryon current in partially twisted, partially quenched chiral perturbation 
theory, we promote Q a from the specific form used in our calculations to the most general 
form transforming under both the adjoint and singlet of SU (7\5). The baryon current has 
the form 




li a D v (B[S„ S u ]BQ a ) + ixpb v (B[S„ S u ]Q a B) + ^D v S U ]B) str(Q a )} 

D„D V (Bv u BQ a ) - D 2 (Bv„BQ a )} + cp \b^D y (Bv v Q a B) - D 2 (Bv v Q a B) 



+ c. 



b^by (Bv u B) - D 2 (Bv u B) str(Q a ) 



(39) 



Restricting all quark indices in Eq. ( |39l) to the sea sector, we can match onto the nucleon 
current of two-flavor chiral perturbation theory in Eq. ( 1381) . Matching with the physical 
light quark charges yields the relations 



C 



gC Q + -C/3 + Cj, 



Cl = -c a 



^Cp, 



(40) 
(41) 



between the partially quenched low-energy constants and the physical parameters of chiral 
perturbation theory. Because our current lacks a flavor singlet component, str(Q a ) = 0, the 
constants /i 7 and c 7 will always be absent from our expressions for nucleon current matrix 
elements. Consequently only the isovector combinations will be expressible in terms of 
physical parameters, specifically fij and cj. The isoscalar combinations will always contain 
unphysical low-energy constants even at the unitary mass point m,- = m u . 



A. Infinite Volume 

A useful check on our formulation of current matrix elements for partially twisted bound- 
ary conditions is the infinite volume limit. In this limit, we must recover the connected parts 
of the proton and neutron form factors. These results, moreover, show the consequences of 
vanishing sea quark charges, and are of use to lattice practitioners beyond the use of twisted 
boundary conditions. 
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FIG. 3: Diagrams contributing to the nucleon vector current in partially quenched chiral perturba- 
tion theory. A thin (thick) line denotes a spin-1/2 (spin-3/2) baryon, while a dashed line denotes 
a meson. Partially quenched hairpins are depicted by a crossed dashed line and the wiggly line 
represents the vector current. 



The calculation of the current matrix elements in Eq. (|36|) can be split into two parts. 
There are local contributions and loop contributions. The local terms are easiest: there 
are Born level charge couplings contained in the free Lagrangian (1131) . and there are addi- 
tional local contributions from higher-order operators appearing in the baryon current (I39|) . 
The loop contributions are generated from the pion-nucleon-nucleon and pion-nucleon-delta 
interactions contained in the Lagrangian Eq. fTLSl) . The relevant diagrams are depicted in 
Figure [31 Additionally at this order, we need to multiply the Born-level couplings by the 
wavefunction renormalization which arises from the diagrams in Figure [TJ 

To express the form factors, we define the three momentum transfer Q = q + B' — B, 
and the quantity 

Here we are additionally considering the nucleon with Fourier momentum transfer q = 
2irn/L, where n is a triplet of integers. Expressing the form factors in terms of Q, we 
can easily generalize between the untwisted case B = B' = 0, and the case of zero lattice 
momentum q = 0. In infinite volume, all results can be expressed as a function of Q. 

For the proton, the local and loop contributions produce the connected part of the electric 
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form factor 



Q2 

G P E (Q 2 ) = 2q u + q d + ^-[ q u (5c a + 2cp) + q d (c a + Ac )] 



6A* 



^jy 2 [dx (2?n + q d ) { \Q 2 log ^ + m%Pf u log 1% } 



(4 

+ f3 ju 



^log^ + m^-SP^logP 2 



(4tt/)2 
+ /?■« 



5 -QHo g y + ml(2-5Pf u )lo g Pf u 



J(m n P n ,A,n) - J(m n ,A,fi) + -x(l - x)Q 2 G(m w P w , A) 



J(m ju P ju ,A,fi) - J(m ju ,A,(i) + - x)Q 2 G(m ju P ju , A) 







L dx - 









(43) 



The coefficients from contributing loop mesons are given by 
for loops containing spin- 1/2 intermediate state baryons, and 



(44) 



Pi = -g(9u-5d), 



= ^ + 2g d ), 



(45) 



for loop containing spin-3/2 intermediate state baryons. From the pion coefficients, one can 
clearly see the photon's coupling to the total charge of the pion, q u — q d . The valence-sea 
meson coefficients, however, reflect that the photon couples to only the valence quarks. The 
connected part of the proton magnetic form factor is given by 



G P M (Q 2 ) = 7: [ Qu(^ a + 2/1/3) + q d {n a + 4/^)] + 



M 



B 



6 



dx 



MbqIn f 1 



47T/ 2 J 

dx [f3' n F(m n P n , A) + (3' ju F(m ju P ju , A)] 



fS^ni^P.^ -\- f3j U 7Tij U Pj U 



(46) 



4ttV 2 Jo 

In writing the above expressions we have made use of abreviations for the non-analytic 
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functions encountered from loop graphs. These functions are 



F(m,5) = —5 log — ^ + VS 2 — m 2 log . , 

V ; & 45 2 B s+VS 2 -m 2 + ie 



m 



G(m,5) = log — 



log 



5 — \Jb 2 — m 2 + i 



is 



45 2 V5 2 - m 2 5 + VS 2 -m 2 Tie' 



(47) 
(48) 



J(m, 5, im) = m 2 log ^ - 2<5 2 log ^ + 25^W^ 2 log * ^ m ' + ' £ , (49) 

and have been renormalized to vanish in the chiral limit. The neutron electric and magnetic 
form factors can be deduced from the above expressions by swapping the electric charges 



G n E (Q 2 ) = G P E (Q 2 ) 



and G n M (Q 2 ) = G P M (Q 2 ) 



Let us focus first on the connected proton form factors at the unitary mass point rrij 



(50) 



Using the physical valence quark charges, we have the connected proton electric form factor 

2 



G p E (Q 2 ) 



2A 2 + 



^ 1 ^|iQ 2 log^ + m 2 P 2 logP 2 } 



(4tt/) 2 

i r 1 

*7) Jo 



9(4tt 



9an 
\4«f) L J<> 



-§Q 2 log^ + m 2 (2-5P 2 )logP 2 
6 n A 



dx 



J(m n P ir , A,/i) - 7(771^, A,fi) + -x(l - x)Q 2 G(m 7T P 7T , A) 



(51) 



Compared to full electric form factor, the connected contribution has the wrong coefficients 
for the tadpole and delta loop contributions. It depends, moreover, on the unphysical low- 
energy constants c a and g 1 , which survive as artifacts of quenching the sea quark charges. 
The situation is similar with respect to the connected contribution to the proton magnetic 
form factor 



G p M (Q 2 ) 



1 

2* 



M. 



B 



36tt/ : 



: (llg 2 A + g A gi + 5£f 2 /4) / dxm 7T P 7T 

Jo 



M B g 2 AN 
24tt 2 / 2 



/ dxF(m n P n , A). 
Jo 

(52) 



Compared to the full magnetic form factor the delta contribution does not have the correct 
numerical factor, and the result depends on unphysical parameters [i a and g±. 

Connected neutron form factors suffer analogous maladies as the reader can easily verify. 
By contrast, the isovector form factors have the correct form. These form factors are defined 
as the difference between proton and neutron form factors 



GUQ 2 ) = G P E (Q 2 )-G E (Q 2 ), 
Gl^Q 2 ) = G P M (Q 2 ) — G n M (Q 2 ) 



(53) 
(54) 
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In the isospin limit, the disconnected operator insertion must cancel out of the isovector 
combinations. Using the connected form factors for the proton and neutron, we find 



G V E (Q 2 



Q 2 

1 + c/-fr + 



O jJL 



+ m 2 P 2 log P 2 

2 ju ju 6 ju 



6(4^ 







L dx< 


9nNN 







+ 9 



juNN 



-§Q 3 log^ + 7^(2-5/*) log J* 
^ 2 log^ + m^(2-5^)log^ 



^9an 
(4tt/) : 



J dxl J{m^P^^ A, //) — J(m w , A, /i) + ^(l — x)Q 2 G(m n P n , A) 



1 

+ 3 



J (m ju P ju , A, //) - J (m ju , A, /i) + -x(l - x)Q 2 G(m ju P jU} A) 



(55) 



where we have abbreviated the combination of couplings 

qInn = ^9a-^9a9i + gf, 

9%nn = 89 a + a 9a9i - gl- 



im) 



These appear as effective axial couplings squared for pion and valence-sea meson loops in 
partially quenched chiral perturbation theory after summing over degenerate mesons. The 
partially quenched isovector magnetic form factor is 



GUQ 2 



Hi 



24vr/ 2 



dx 



glNN m nPn + g 2 juNN m juPju 



M B g 2 AN 

12n 2 f 2 



dx 



F{m«P«, A) + -F{m ju P ju , A) 



(57) 



These partially quenched form factors agree with those determined in SU(4\2) partially 
quenched chiral perturbation theory [2?], 36]. The local contributions are now proportional 
to fii and c/, which are physical parameters. Both of these partially quenched form factors, 
however, depend on the unphysical coupling g\. Taking the valence-sea meson to be degen- 
erate with the pion, m 2 u = m 2 ^ , this dependence disappears because g 2 NN + g 2 U NN = ^9 ji- 
lt is only in this limit that the isovector form factors reproduce the correct QCD physics. 3 



3 This point is often overlooked, particularly in mixed action simulations which are automatically partially 
quenched. In a mixed action simulation, the valence and sea pion masses are tuned in order to mitigate 
unitarity violations. The valence-sea meson mass, however, is not protected from additive renormalization 
and is degenerate with the pion only in the strict continuum limit. 
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FIG. 4: Additional diagrams contributing to the nucleon vector current in finite volume PQxPT. 
Diagram elements are the same as those in Figure [3j In the infinite volume limit, these diagrams 
vanish by Lorentz invariance. 



For completeness, the nucleon isovector form factors resulting from taking mj u = m 2 are 

GUQ 2 ) 



Q 2 2 

1 + Cl-TT + 



K (^f) 2 Jo 



dx 



^log^ + m^logP 2 



2g 



(4tt/) 2 
%9an 



dx 



~Q 2 log^ + m 2 (2-5P 2 )logP 2 



dx 



3(47T/)2 

for the isovector electric, and 



J (m n P n , A, fx) - J (77V, A, n) + -x(l - x)Q 2 G(m n P n , A) 



G V M {Q 2 



9 2 a Mb 

27T/ 2 J 



I dxm^P^ - ~^§jf [ dxF(m n P n ,A), 
Jo y7r / ^0 



(58) 



(59) 



for the isovector magnetic form factor. These results agree with the standard two-flavor 



chiral perturbation theory calculations in the literature [37|, [38 



B. Finite Volume 

We now evaluate the matrix elements contributing to the connected part of the proton 
current in Eq. (|36|) in finite volume. This requires us to revisit the computation of the 
wavefunction renormalization diagrams shown in Figure HJ and the form factor diagrams 
shown in Figure [31 Additionally there are new contributing diagrams which are displayed 
in Figure HI These diagrams ordinarily vanish in infinite volume by Lorentz invariance. 
Furthermore at finite volume with periodic boundary conditions, these diagrams also vanish 
but by the remnant discrete rotational symmetry (cubic invariance). With continuous twist 
angles, however, these diagrams do not vanish and are required in our computation of current 
matrix elements. 

Resulting expressions for the temporal and spatial components of the current are quite 
lengthy and are displayed in their entirety in Appendix El For ease, the expressions given in 
this section will employ various simplifications. Firstly we work at the unitary mass point, 
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m ju = m n- We will focus on the connected proton result, as well as the isovector combination 
of finite volume matrix elements. Additionally as Lorentz symmetry is not respected at 
finite volume, the form factor decomposition in infinite volume is no longer valid, see [39], 
for example. With twisted boundary conditions, we find the temporal component of the 
current acquires spin dependence at finite volume. Similarly the spatial components of the 
current acquire spin diagonal terms. These terms are displayed in Appendix [AJ while the 
expressions presented here will either be unpolarized for the temporal component, or the 
difference of polarized matrix elements in the case of the spatial components. Lastly the 
results in Appendix |A] are for a general frame of reference in which the initial state moves 
with momentum 8/L, and the final state moves with 6'/L. The expressions given in this 
section will be specific to either the rest frame, in which 6 = 0, or the Breit frame, in which 
= -0'. 



1. Rest Frame 



In the rest frame, the momentum transfer is given by Q = q + B' . The finite volume 
modifications to proton current matrix elements are given by 



1 x 1 f l 

-^2(p,m\5J 4 \p,m) = — dx 

1 m=± J Jo 



I 1/2 {m 7T P 7T , xQ) - ^Li/ 2 (m n , 0) - ^Ti^m^, B') 



12f 



\\g\ + g A9l + 5^/4) J(m n , 0, 0) + J(m„, B\ 0) 



30ajv [J(m^ 0, A) + J(m^ B' , A)] 



1 



dx 



(llg 2 A + g A9l + hgl/A)J{m^ 0, Q, xQ, 0) 



6/ 2 Jo 

0,Q,xQ,A) 



(60) 



for the unpolarized time component of the current; and, 



(p,±\5Ji\p,T) 



-i(±| [s k ,Sj] |t) 1^(13^ + 11^1 + </?/4)x? (7^,^,0) 



5 ki ^ L K, j {m 7n B',A) 



6 



+ 7:Qfc / dx 
6 Jo 



(llg 2 A + g A g x + 5gl/4)&\m n P n , 0, Q, xQ, 0) 



+ -gl N a\m n P 7T ,0,Q,xQ,A) 



(61) 



for the spatial components. We have chosen spin-flip matrix elements; these are simply 
related to differences of spin polarized matrix elements. 
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The finite volume corrections to isovector matrix elements, we write out similarly. 



-\2(p,m\5Jl\n,m) = — / dx 



2/2 



9l 



l 1/2 (m n P n ,xQ) - -^(m,, 0) - ^X 1/2 (m n , B') 



J(m w ,0,0) + J(m w ,B',0) 



9 



9an 

[J(m n , 0, A) + J(m n , B', A)] 



g 2 A J(m n P n , 0, Q, xQ, 0) - -g\ N J{m^ 0, Q, xQ, A) 



(62) 



for the unpolarized time component of the current; and, 



(p,±\5J+\n, T ) = 7 i(±\[S k ,S J ]\T){25 kt (g 2 A + g A g 1 )W(m 7T ,B',0) 



P 



+3Q k / dx 
'o 



g 2 A V\m„P n , 0, Q, xQ, 0) + -g 2 AN &\m n P n , 0, Q, xQ, A) 



(63) 



for the spin-flip spatial current. 



2. Brett Frame 

In the Breit frame, we choose B' = —B and the the momentum transfer is thus given by 
Q = q — 2B. The finite volume modifications to proton current matrix elements are given 
by 

i x i r 1 

- (p, m\5J±\p, m) = — \ dx \T l/2 (m^, xQ + B) - ZWm*, B)] 

m=± f J 

{llg 2 A + g A9l + 5^/4) J(m n , B, 0) - 3g 2 AN J{m n , B, A) 



6/ 2 



6/ 2 Jo 



dx 



(llg 2 A + g A9l + 5g 2 /4)J(m n P n , B,Q + B,xQ + B, 0) 



B,Q + B,xQ + B,A) 



(64) 
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for the unpolarized time component of the current; and 
1 „, , . f .2 

7 



(p,±\6Ji\ Pl T) = -^(±\[S k ,S,]\T)i -5 k ^(l3g 2 A + llg A g 1 + g 2 1 /4)W(m n ,B,0) 



9 

+iQk I dx (llg 2 A + g A9l + 5g 2 1 /A)a i (m 7T P 7T , B,Q + B,xQ + B } 0) 
6 Jo 



+ -g 2 AN JC ji (m n P m B,Q + B,xQ + B, A) 



(65) 



for the spatial components. 

The isovector current matrix elements can similarly be derived in the Breit frame. For 
the time component of the current, we have 

1 x 1 r 1 

-y2(p,m\5J£\n,m) = — / dx [l 1/2 {m n P n , xQ + B) - X 1/2 (m 7r , B)] 

Z m=± J J 



3 

T 2 

/ 2 



(^(m*, B, 0) - -^J(m^ B, A) 



+— / <ir 
'o 



g 2 A J(m n P n , B,Q + B,xQ + B, 0) 
-^g 2 AN J{m^,B,Q + B,xQ + B,A) . (66) 



Finally, the spatial isovector current has the spin-flip finite volume corrections given by 

1 



(p,±\5Jt\n, T ) = -^(±1 [£*,£,•] |=F>^ - 4fe(^i + ^i)/C J (m^,B,0) 



+3Q fc / dx g 2 A C ji ( m7T P w ,B,Q + B,xQ + B,0) 
Jo L 

+ \gl N C ji {m^ B,Q + B,xQ + B, A) 



(67) 



C. Numerical Estimates 

To estimate the effect of finite volume corrections, we use phenomenological input for 
the various coupling constants. The values we use have been listed above in Section IHIl 
We restrict our attention to isovector quantities and numerically evaluate the corrections in 
the rest frame. We will comment on the qualitative behavior of volume corrections in Breit 
frame. 

Consider first the finite volume corrections to the isovector electric form factor G V E (Q 2 ). A 



4 Strictly speaking there are no longer electric and magnetic form factors on a torus as the decomposition 
in Eq. ([30]) relies on Lorentz invariance. We will use electric (magnetic) to denote quantities calculated 
from the temporal (spatial) component of the current with the appropriate spin structure. 
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FIG. 5: Relative change in the isovector electric form factor due to twisted boundary conditions 
in the rest frame. Plotted versus the twisting angle & is AG V E (Q 2 , L) given in Eq. (|69p . The lattice 
size L is fixed at 2.75 fm, and the momentum transfer is |Q| = O'/L = 6' x 0.072 GeV. 



To access this form factor, we use unpolarized matrix elements of the time component of 
the current. We find 

GUQ 2 , L) = GUQ 2 ) + h[GUQ% (68) 

where G V E (Q 2 ) is the infinite volume form factor given by Eq. fl58|) . and 5l[G v e (Q 2 )] is the 
finite volume correction which is identical to the unpolarized matrix element in Eq. (1621) . 
Here we work with the momentum transfer entirely due to twisting Q — B' — 6'/L, and 
take 0' to lie along one spatial direction. Notice the finite volume correction to the isovector 
electric form factor is independent of any unphysical parameters, in particular the coupling 
g\. The infinite volume isovector electric form factor depends on the parameter c/(/x), the 
value of which can be inferred from the charge radii of the proton and neutron. Using the 



Particle Data Group averages [40], we find c/(/i = 1 GeV) = —0.393. In Figure [5j we plot the 



relative change in the isovector electric form factor due to volume effects AG V E denned by 

AG E (Q , L) . (69) 

Here we keep the box size fixed at 2.75 f m, and plot versus the twisting angle 9'. Qualitatively 
the finite volume effect oscillates about the infinite volume form factor as 9' is increased. 
The oscillations are damped, but this behavior is apparent at momentum transfers too large 
to trust the effective theory. There are no finite size effects at 9' = because of charge 
non-renormalization (which holds due to treating the time direction as infinite). The results 
are shown for = 0.25 GeV with the finite volume effect generally at the percent level or 
less. The effect of finite volume is of course smaller for larger pion masses. 

Considering the spatial components of the current, we can determine the magnetic form 
factor. Additionally there are volume corrections to the spatial current, and the net effect 
has the form 

G\ d {Q\ L) = G V M (Q 2 ) + S L [GUQ% (70) 

where G V M (Q 2 ) is the infinite volume form factor given by Eq. (159]) and 8l[G v m (Q 2 )] is the 
finite volume correction, which follows from Eq. (1631) . Choosing for simplicity B' = B'y, 
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FIG. 6: Relative change in the isovector magnetic form factor due to twisted boundary conditions 
in the rest frame. Plotted versus the twisting angle 9' is AG V M (Q 2 ,L) given in Eq. (|72p . The lattice 
size L is fixed at 2.75 fm, and the momentum transfer is |Q| = 9'/L = 9' x 0.072 GeV. The bands 
arise from uncertainty in the low-energy constant g%, which we vary assuming it is of natural size. 



and utilizing the z component of the current between an initial state spin-up and final state 
spin-down, we have 



Sl[Gm{Q 2 )] 



g M p + 9Agi)JC 2 (m n , B'y, 0) 
3M N 



P 



dx 



g 2 A C 33 (m w P n ,0, B'y,xB'y,0) 



+ yg 2 AN £ 33 (mnPn, 0, B'y, xB'y, A) 



(71) 



Notice this volume correction depends on the unphysical coupling gi which arises as a 
consequence of having enlarged the valence flavor group. The infinite volume isovector 
magnetic form factor depends upon the parameter /// which we can estimate using the known 
values of the proton and neutron magnetic moments. We find [ij = 6.77. In Figure [61 we 
plot the relative change in the isovector magnetic form factor due to volume effects AG V M 
defined by 

G V M (Q 2 ) 



G v m(Q 2 ,L) 



(72) 



23 



m n = 0.25 [GeV] 
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FIG. 7: Comparison of finite volume and infinite volume isovector magnetic form factors. Plotted 
as functions of Q 2 are the infinite volume form factor G V M (Q 2 ), and finite volume form factor 
G V M (Q 2 ,L). The lattice size is 2.75 f m, and the value of the unknown axial coupling has been fixed 
to gi = —0.5. 



Again we keep the box size fixed at 2.75 fm, and plot versus the twisting angle 6'. Because 
the effect is non-negligible, we choose a few values of the pion mass. The result, moreover, 
is sensitive to the value of gi which has been varied assuming natural size, —2 < g\ < 2. In 
Figured we compare the extracted form factor at finite volume G V M (Q 2 , L) with the infinite 
volume form factor G V M (Q 2 ) as a function of Q 2 = 6' 2 /L 2 . In this figure, we keep the lattice 
size at 2.75 f m, and fix the pion mass to be 0.25 GeV. Furthermore, we choose the value of gi 
favored by comparing with SU(3) chiral perturbation theory, namely g\ = 2(F—D) « —0.5. 5 
Lastly we comment on the size of volume corrections in the Breit frame. Comparing 
the expressions for the isovector electric form factor in the rest frame Eq. (162]) . and the 
Breit frame Eq. flrJBl . we see that all factors depending on B but not Q are doubled in the 
Breit frame. This is due to the symmetry under the exchange of the initial and final state 
twists: the finite volume functions are even. Because there is some cancellation among the 
contributions to the finite volume electric form factor in the rest frame, we can anticipate 
that the finite volume corrections in the Breit frame will generally be of the same size. For 
the magnetic form factor, comparing Eq. (I63I) and Eq. (|67|) shows similarly that the effect 
from the K? terms doubles. While this function is odd with respect to argument, terms 
from the initial and final states add coherently because there is a relative sign from the spin 
algebra. Because empirically we observe the dominant volume correction arises from the K? 
term, the volume effect for the magnetic form factor will roughly double in magnitude in 
the Breit frame. Given that the coefficient of this term depends upon the unphysical and 
unknown parameter gi, the Breit frame does not offer an advantage over the rest frame. 
A lattice calculation of g\ is necessary to control the systematic uncertainty from volume 



effects in this approach. This is not the case for isospin twisted boundary conditions |12 
see Appendix 03] 



One could calculate g\ directly by determining the axial couplings of hyperons in the SU(3) limit (and in 



the chiral regime). The first lattice calculation of hyperon axial charges has been recently performed 41 1 
but naturally with a focus on SU (3) breaking. 
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V. SUMMARY 



In this work, we compute finite volume modifications induced by partially twisted bound- 
ary conditions. We utilize heavy baryon chiral perturbation theory in finite volume. Baryons 
are embedded into representations of S77(7|5), where the extra flavors are fictitious, and dif- 
fer only in their boundary conditions. The nucleon mass splittings are determined, and 
demonstrated to be negligible on current-sized lattices. The main focus of our work is the 
derivation of finite volume corrections to the vector current matrix elements of the nucleon. 
Continuous momentum is inserted on the active valence quark lines using flavor changing 
currents in the enlarged flavor group. Disconnected operator insertions cannot be accessed 
at continuous momentum using this technique, and our calculation is therefore restricted to 
connected current insertions. Isospin breaking and cubic symmetry breaking lead to vari- 
ous structures not encountered in infinite volume. We give complete expressions for finite 
volume current matrix elements using general kinematics. To estimate the size of these 
corrections, we choose rest frame kinematics, and consider both the spatial and temporal 
components of the current. Generally the volume corrections lead to oscillatory behavior 
about the infinite volume answer. In the region of small twist angles, the volume effects 
can become rather pronounced due to terms that break cubic symmetry. To extract the 
isovector magnetic moment and electromagnetic radii from lattice data at zero Fourier mo- 
mentum, careful determination of volume effects will be required. This is complicated by 
the dependence on an unphysical and unknown axial coupling g\. As shown in Appendix [Bj 
a different implementation of twisted boundary conditions can eliminate this dependence. 
In this implementation, there are no fictitious flavors, rather the isospin transition is simu- 
lated directly. Compared to the meson sector, the baryon sector appears more susceptible 
to volume corrections due to partial twisting. The Breit frame kinematics do not simplify or 
reduce the volume corrections. Partial twisting provides a novel way to probe any isovector 
nucleon matrix element at continuous values of momentum transfer. The formalism de- 
veloped here allows one to compute finite volume modifications to these observables, and 
thereby control the extraction of moments, radii, etc, from lattice QCD data. 
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APPENDIX A: FINITE VOLUME CURRENT MATRIX ELEMENTS 

In this Appendix, we list the finite volume corrections to current matrix elements. For 
ease of presentation, we remove the Pauli spinors. Here we work in a general frame where the 
initial-state nucleon has momentum B = 6/L, and the final-state nucleon has momentum 
B' + q, where q is a Fourier momentum mode of the lattice and B' = 6'/ L. Here we list 
only the proton matrix elements as a function of q u and q^. One can use charge symmetry, 
qu q<ii to deduce the neutron matrix elements. As explained in the main text, each result 
includes only connected part of the matrix elements. 
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The finite volume modification to the time-component of the current matrix element in 
Eq. fl36|) reads 



<5J 4 



■ Id) 



3 

172 



Qd + ~Qu 



T 1 / 2 (mj u P ]u ,xQ + B) - -Z x / 2 (mj U ,B) - ]-l 1/2 (m JU , B') 



9"iY 3 iV3^(mvr,0,0) +g] uN3N3 J(m ju ,O,0) + -g? NaNa (J(m*, B, 0) + J(m*, B' , Q)J 



-<-iu 



9'juN a N 3 (J{m ju ,B, 0) + J(m ju , B', 0) 



s'iViiV^C"^, 0, 0) + 9 2 juNlNl J(m ju ,0, 0) + -g'J NlNl (j(m^,B, 0) + J(m^, B' , 0)) 



jS^NiiV! {J(mju,B, 0) + J(m JU , B', 0)) 



5 /— , 



J(m w ,0, A) + 2J(m ju , 0, A) + - ^(m^, B, A) + ,7(m T , B', A) 



+ 2 ( J(m ju , B, A) + J{m ju , B' , A) J 



2/ 2 



A) + 2J r (m i „, 0, A) + - (J(m^, B, A) + J(m w , B', A) ) 



^iV 3 iV3^( m -. 0. 0) + gtel - 9A91 - gi/2) [J(m w , B, 0) + J(m n , B' , 0)J + gj uNaNa J{m ju , 0, 0) 



qI^N^^, 0, 0) + g (Sa + %A9i - S?/2) (JK, B, 0) + J(m^, B', 0)J + <?f uJVlJVl J(m ju ,0, 0) 
<7d + J ( ^J(m^, 0, A) + i (j(-m^, B, A) + ^(m^ B', A)) + ^J(m ju ,0, A) 



+ §g u (ij(mi,0, A) + -(j( miT ,B, A) + J(mn,B', A)) + ^(m ju ,0, A) 

~[S-Q,Sj] f dx[f3 7r J ] (m 7T P 7r ,B,xQ + B,0) + l3 ju J : >( mju P JU ,B,xQ + B,0)] 

I Jo 

_!^v [ S . q, 5j ] y dx [/3' n ji(m^P^, B, xQ + B, A) + 0' ju J J (m ju P ju , B, xQ + B, A)] 

-JL /" dx [Ar .7(711^ P T ,B,Q + B,xQ+B,0)+/% u J(m JU P jtl ,B,Q + B,:rQ + B,0)] 
*J Jo 

3<?L 



- / dx[p!rJ(m n P n ,B,Q + B,xQ + B,A)+P' ju J(m jn Pju,B,Q + B,xQ + B,A)]. 
Jo 

The effective axial couplings have been given above in Eqs. fl2Tj) and fl23|) . while the loop 
coefficients appear in Eqs. (1440 and f|45|) . The spatial components of the current matrix 



(Al) 
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element in Eq. ( 136]) receive the finite volume modification 

SJi = -j^l («d + ifc.) 



+ |«u [9?N lNl B, 0) + K*(in*, B',0)) + flj^^ (?C l (m Jtl , B, 0) + (wy ttJ B', 0))] J 



[Si , (Sj 



9d + \9'Jn 3 N3 • B ' °) " Ki ("** ' B '> °) ) + ^'iVaJVa K«. B, 0) - /C J (m JU ,B',0) 



+ \lu [g? NlNl (K?{m„, B, 0) - K^m^B', 0)) + g£ NlNl (k? (rrij u , B, 0) - tf(m itlJ B', 0))] J 
- 9u J [- (rK, B, A) + ^(m^, B', A)) + - (lC{m ju ,B, A) + K^m^, B', A)) 



-3^1 I"'' 2 



|g„J ^( A ' .-B.Ai- /v'' ( /(/•..- . B'. A i 



+9AAT 



3/2 



5d - 



i?u j [J (iCi (m n , B, A) - K. ] (m„ B', A)) + | (/C J (m ju , B, A) - /C J (m JIt , B', A) 



"472 / dx [^C i (m w P w ,B,Q + B,xQ + B,0) + /3 ju C i (m ju P ju ,B,Q + B,xQ + B,0)] 

+ ^rjf t dx [AC l (m^P^, B,Q + B,xQ+ B, A) + 0. u C (m 3 u P 3 , u , B,Q + B,xQ + B, A)] 
Jo 

[Q ■ S, Sj] f dx [p 7C C ji (m v P vl B,Q + B, xQ + B, 0) + /3 :ju C : > i (m ju P :ju , B,Q + B,xQ + B, 0)] 
Jo 

[Q S, Sj] [ l dx [/^^(m^ , B, Q + B, xQ + B, A) + p' ju C : > i {m jv .Pj u , B, Q+B,xQ + B, A)] 
Jo 



3 
2/2 

2/2 



(A2) 



Appearing in the above expressions for finite volume modifications are functions depend- 
ing on the difference of finite volume mode sums and infinite volume momentum integrals. 
The various definitions are as follows: 



l 1/2 (m,A) = jaX] 



dk 



J{m,A,B,C,A) = dXX 



L 3 „ [{k + A) 2 +m 2 ]V2 J ( 27r )3 [ fe 2 + M 2]l/2' 

dk {k + A)-(k + B) 



1 ^ (k + A)-(k + B) 

Ifi 1^ [(fc + C)2 + 02]6/2 " j (2^) 3 [(fe + C) 2 + /3|] 5 / 2 



J J (m, A, B, A) = / dAA 
'0 



r,3 



(fc + A)j 



(fe + A) 



/■OO 1 

tf(m,B,A) = ^ dA IsE 
£ J '(m,A,.B,C , ,A) = 



£ 3 V [(fe + B) 2 +/3|] 5 / 2 7 (2^) 3 [(fc + B) 2 +/?|] 5 / 2 

(fe + By 



[(fc + B) 2 + /3l] 3 / 2 ' 



(A3) 
(A4) 

(A5) 
(A6) 



dX 



1 ^{k + A) ■ (k + B)(2k + A + B)' J 



I? 



[(fc + C) 2 + /?|] 5 / 2 
dk (k + A)-(k + B)(2k + A + By 



2tt) 



[(fe + C) 2 + /3lp/2 



(A7) 



and 



C J {m,A,B,C,A) = / dX 
Jo 



1 x - (fe + A) l (2k + A + By 
[(fe + C) 2 +/3|] 5 / 2 



dfc (fe + A) l (2fc + A + B) 
(2tt)3 [(fc + C) 2 + /3i] 5 / 2 



(A8) 
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We also use the short hand J{m, A, A) = J(m, A, A, A, A). We show how to evaluate 
these functions numerically in Appendix O 



APPENDIX B: FINITE VOLUME ISOVECTOR CURRENT MATRIX ELE- 
MENTS FROM 5J7(4|2) 

In this Appendix, we detail the finite volume corrections to current matrix elements using 
the alternate implementation of partially twisted boundary conditions proposed in 0, 12 



This method does not rely on fictitious flavors of valence quarks. Instead, one directly 
confronts the isospin changing operators whose matrix elements give rise to isovector form 
factors. From the outset, one is aware that disconnected contributions cannot be accessed. 
The flavor structure, moreover, only requires a simple modification of existing partially 
quenched theories to include twisted boundary conditions. Results for the finite volume 
isovector magnetic form factor under simplifying kinematics were given in however, as 
complete expressions for finite volume current matrix elements were not presented, we give 
the complete expressions here. 

Let us briefly summarize the setup used in 12{|. We restrict our attention to an S77(4|2) 



theory with quarks contained in a field Q, which is given by Q = (u,d,j,l,u,d) T . Each 
quark is periodic but coupled to a uniform Abelian gauge potential of the form = 
diag( J B^, J B^,0,0, J B;, J B^), with B^ = (6 U /L,0) and B d = (0 d /L,O). In this formulation, 
momentum is injected by isospin changing operators provided B d 7^ B d . Keeping these 
twists different introduces isospin breaking via finite volume effects. The partially quenched 
isospin splittings for the pion were numerically demonstrated to be quite small on current 
lattices [l3|. To calculate the nucleon isospin splitting, we evaluate the sunset diagrams 
shown in Figured] in the partially twisted S77(4|2) theory. The nucleon isospin splitting is 
given by 

M n -M p = -^l~g] uNN [K:{m ju ,B d ,0) - K,(m ju , B u ,0 

- \qIn [JC(rn ju , B d , A) - K(m ju , B u , A)] J . (Bl) 

The effective axial coupling gj uNN has been given in Eq. (1561) . When the twists are isospin 
symmetric, the nucleon mass splitting accordingly vanishes. The maximal isospin splitting 
is occurs when B = 1, 1) for one flavor, and B = for the other. On current lattices 
this maximal splitting is at the percent level and can practically be ignored. 

For the operator J+ = u'~f f j - d, continuous three-momentum of the form B n = B u — B d is 
induced in flavor changing matrix elements. Thus we consider the isovector-vector current 
matrix elements between nucleons 

(p(q)\4H°)) ^ (P(P')\J>(P)) = (P(P')\J7\P( P )) - (n(P')\j; m HP)}. (B2) 

On the left, we have denoted only the Fourier momentum. On the right, the momentum of 
the initial-state nucleon due to twisting is P = B u + 2B d , while the final-state nucleon has 
momentum P' = q + 2B U + B d . The momentum transfer we denote by Q and is given here 
by Q = q + B n . The equality between the isovector-vector current and differences of the 
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electromagnetic current matrix elements follows from the SU(2) va \ ence symmetry subgroup 
of the full S77(4|2) group. At finite volume, this symmetry is broken and one must address 
the volume corrections to the matrix element on the right-hand side. To determine these 
corrections, we evaluate the one-loop diagrams for the isospin transition matrix element 
using the partially twisted SU(4\2) theory. The relevant diagrams are shown in Figures [3] 
andU 

The finite volume modification to the time-component of the isovector current matrix 
element in Eq. (IB2j) reads 



SJ4 = / dx 
J Jo 

8/ 2 ' 



Zi/2(mjuPju, xQ + B d ) - ^Ii/2(m.j u ,B u ) - ^l 1/2 (m ju , B d ) 



jNN 



J(m^, 0,0) + J(m^,B n ,0) 



2 

■ 9juNN 



J{m ju ,B u ,0) + J(m ju ,B d ,0) 



9 AN 



6p 



3J(m„,0, A) + 3J(m w ,B w , A) + J(m ju , B u , A) + J(m ju , B d , A) 



1 r 1 

— [Q ■ S, Sj] / dx [gl NN J 3 KP„0, xQ, 0) + g] uNN J J (m JU P JU , B d , xQ + B d , 0) 
J Jo 

[Q ■ S, Sj] f dx 
Jo 



9 an 



J 3 KP„0, xQ, A) + -J 3 (m JU P JU7 B d , xQ + B d , A) 



1 r 1 

— / dx [g* NN J(m n P n ,O,Q,xQ,0) + g] uNN J(m ju P ju ,B d ,Q + B d ,xQ + B d ,0)] 
J Jo 



4/ 



9 a n 



P 



dx 



J(m n P^,0,Q,xQ,A) + -J(m ju P ju ,B d ,Q + B d ,xQ + B d ,A) 



(B3) 



We have omitted writing the Pauli spinors here and below. The effective axial couplings, 
g^NN an d g'juNNi have been given above, Eq. fl56l) . The spatial components of the current 
matrix element in Eq. fl36l) receive the finite volume modification 



SJ+ = 



1 1 



2/ 2 [6 

[Si , Sj 
~6p~ 



~^9juNN 



K} (m 3U ,B u ,0) + JC l (m ju ,B d ,0) 



2 2 

g9AN 



IC l (m JU ,B u ,A)+IC l {m JU ,B d ,A) 



gl NN K. 3 (m„,B„,0) + g? NN IC J (m ju ,B u ,0) - K? (m ju , B d , 0) 



+ 



9 AN [Si> Sj 



3/ 2 



K?(m*,B*,A) + - [V(m ]u , B u , A) - V {m ]u , B d , A) 



1 f 1 

I dx [gl NN C{m^P 7 ,,0,Q,xQ,0) + g] uNN C(m iu P lu ,B d ,Q + B d ,xQ + B d ,0)] 
J Jo 



8/ 



9 AN 



2/ 2 Jo 



dx 



C(m n P w ,0, Q, xQ, A) + -jC\m ju P ju ,B d , Q + B d , xQ + B d , A) 



1 f 1 

-— 2 [Q • S, Sj] / dx [g* NIf £? % (m«P K ,0, Q, xQ, 0) + g 2 3uNN C : > l {m JU P JU ,B d , Q + B d ,xQ + B d , 0) 
4 J Jo 



9 AN 

2p 



[Q ■ S, S 3 ] f dx 
Jo 



C ji (m w P n ,0, Q, xQ, A) + -C ji (m ju P ju ,B d , Q + B d , xQ + B d , A) 



(B4) 



From these expressions, we can simplify things by forming unpolarized (polarized) ma- 
trix elements for the temporal (spatial) part of the current. Furthermore, we restrict our 
attention to the unitary mass point, where m? = m%, and choose the twist parameters such 
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that 6 d = and O u = 6, which corresponds to rest frame kinematics. 6 The time component 
of the current becomes 



- M<w 4 + M 



m— ± 



P Jo 
3 



dx 



9 2 A 



l 1/2 {m 7r P n ,xQ) - -li/ 2 {m w ,0) - ^X 1/2 {m w ,B) 



J(m- K ,O,0) + J(m^,B,0) 



1 2 

^9an 



J{m„, 0, A) + Jlrn^B, A) 



"T 5 ' 



J(m„P„, 0, Q, zQ, 0) - -g 2 AN J(m„P„,0, Q, xQ, A) 



(B5) 



with Q = q + B. The spatial current reads 



1 



±WW) = -7o(±\[Sk,Sj}\ T ){26 ki 



P 
+3Q fc 



g\ K? (m» , B, 0) + A) 



^ ^{m^, 0, Q, xQ, 0) + -g 2 AN O l {m„P^ 0, Q, xQ, A) 

y 



(B6) 



We have chosen spin flip matrix elements; these are simply related to differences of spin 
polarized matrix elements. Because the the finite volume modifications proportional to K? 
are non- vanishing only in the directions with non- vanishing twist, the spin and momentum 
transfer structure of these terms are identical to the magnetic part of the current matrix 
element. Consequently one cannot be sensitive to the magnetic form factor without addi- 
tionally acquiring finite volume modifications from K? terms. These terms are seen to be 
numerically larger than £ J \ especially for small twists Finally, notice these results are 
independent of the unphysical parameter g\. 



APPENDIX C: FINITE VOLUME SUMS 

In this Appendix we describe the evaluation of the mode sums required for the finite 
volume corrections to the nucleon mass and isovector form factors above. In the main text, 
we have used various functions entering in the computation of loop graphs in finite volume. 
Here we evaluate each function systematically in terms of Jacobi elliptic functions and error 
functions. 

The basic sums required are of the form 

p{Q,M) - -jj 2_j ^ _j_ + _A/f2j/3 ~ J (2tt) 3 [k 2 + M 2 f 

ji , , i v k { r dk w 

p{q,M) ^Lp + ^+^Y J (27r) 3 [(fc + q) 2 + A4 2 ]/ 3 ' 

} ~ L 3 ^ [(k + q) 2 + M 2 ]P J {2-nf[{k + q) 2 + M 2 Y' { } 



6 The choice 8 d = —9 U does not result in any dramatic simplifications. In particular there will be residual 
gi dependence in this case. 
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The latter two functions can be derived from the first via differentiation, explicitly the 
relations are 

Tp(q,M) = ^ £ i Tf l - 1 (q t M) - q%(q,M), (C2) 



and 



n{qM) 



2{fi - 1) dq l 

1 d 2 
4(f3-l)(f3-2)d q ^- 2{q ' M) 



+ 



S ij + q l ^~ + q°^)j lp-i(q,M) + q i q%{q,M)- (C3) 



203-1) 

Evaluating the first function in Eq. (jCip for arbitrary (3, we find 

I poo 3 

■lAq.M) 



87r 3 / 2 r(/3) 



dTT 0-* e -rM' 



£1 

4t 



(C4) 



where ^3(9, 2) is a Jacobi elliptic function. 

In the main text, we utilized several different mode sums in the evaluation of finite volume 
effects. We now write them out in terms of the basic finite volume functions Xp(q,A4), 
Z l p{q, M), and lp(q, M). With the abbreviation p\ = m 2 + 2AA + A 2 , specifically we have 



J(m, A, B, C, A) 

J j (m, A,B,A) 
fC(m, A, A) 
/C J '(m, A, A) 
C ij (m, A, B, C, A) 



dAA 



^ /2 (C,/3 A ) + (A + B)% /2 (C,p A ) + A • £Z 5/2 (C,/3 A ) 



-5/2 



dXX 
dX 
dX 
dX 



li /2 (B,P±) + A% /2 (B,p A ) 



8"T' /2 {A, p A ) + 2A*T 3/2 (A, p A ) + A 2 1 3/2 (A, p A 

li /2 (A,p A ) + A% /2 (A,p A ) 

2T 5 j /2 (C,p A )+2A% 2 (C,p A ) 
+ (A + Byr 5/2 (C,p A ) + A\A + Byi 5/2 (C,(3 A ) 



(C5) 
(C6) 

(C7) 
(C8) 



(C9) 



Lastly the evaluation of the A-integrals can be done in closed form. For completeness the 
required A-parameter integrals are 



(1Xr -r(X> + 2 b x + <?) = I /V(* 2 - C 2 )Erfc(&v^), 
2 V t 



dX X e -^ 2 + 2b ^) 



2V r 



- 6e r62 Erfc(&^) 



(CIO) 
(Cll) 



where Erfc(x) = 1 — Erf(x), and Erf(x) is the standard error function. 
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